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basic  research  in  missile  control  theory  for  potential  system 
applications . 
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SECTION  I 


INTRODUCTION 


1.  OBJECTIVE 

This  report  documents  an  autopilot  design  for  a  tactical 
bank-to-tum  air-to-air  missile  that  must  intercept  a  highly 
maneuverable  target.  The  approach  used  in  this  study  takes  into  account 
the  aerodynamic  and  inertial  interactions  (cross  coupling  effects)  that 
are  common  to  this  type  of  problem.  This  approach  is  different  from 
many  approaches  to  air-to-air  missile  autopilot  design  since  some 
effects  of  coupling  (primarily  inertial)  are  often  neglected  during  the 
design  stages  of  the  autopilots.  Such  neglecting  of  terms  leads  to  a 
simplified  nonlinear  model.  The  approach  taken  for  the  project  is  to 
develop  the  linearized  model  from  a  full  nonlinear  model.  A  stochastic 
approach  is  applied  to  gain  scheduling  which  allows  large  random 
parametric  variations,  this  approach  leads  to  simplified  gain  scheduling 
requirements.  It  is  expected  that  this  approach  will  enhance  the 
performance  of  the  autopilot  in  the  presence  of  sensor  noise  and 
modeling  errors.  Improved  performance  leads  to  decreased  miss  distances 
during  the  missile- target  intercepts. 

2.  TERMINOLOGY  AND  REFERENCE  FRAMES  USED 

Figure  1  shows  the  missile  body  coordinate  frame  which  is  used  to 
describe  the  six-degrees-of-freedom  (6  DOF)  for  a  missile.  Three- 
degrees-of-freedom  are  associated  with  the  missile  body  position  and  the 
other  three-degrees-of-freedom  are  associated  with  the  missile  body 
angular  orierrat ion. 

The  missile  body  center  of  mass  is  the  origin  for  the  coordinate 
frame.  From  Figure  1  it  should  be  apparent  that  the  x  axis  passes 
through  the  center  of  mass  of  the  missile  from  the  rear  to  the  front. 

The  y  axis  passes  through  the  certer  of  mass  of  the  missile  from  the 
left  side  of  the  missile  to  the  right  side  of  the  missile.  The  z  axis 
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Figure  1.  Missile  Body  Coordinate  Frame 


passes  through  the  center  of  mass  of  the  missile  from  the  top  of  the 
missile  to  the  bottom. 


The  missile  body  velocity  components  u,  v,  and  w  lie  along  the  x, 
y,  and  z  axis,  correspondingly.  The  right  hand  rule  defines  the  sense 
of  the  rotation  about  the  three  axes .  The  angular  rotations  about  the 
x,  y,  and  z  axes  are  correspondingly  called  roll,  pitch,  and  yaw  angles. 
The  term  body  rates  is  used  to  describe  the  rotation  rates  (roll,  yaw, 
and  pitch  rates)  about  the  x,  y,  and  z  axis. 

Three  planes  can  be  defined  from  Figure  1.  The  y  and  z  axis  form 
the  roll  plane.  The  x  and  y  axis  form  the  yaw  plane.  The  x  and  z  axis 
form  the  pitch  plane.  The  angle-of-attack  (a)  is  the  angle  formed  in 
the  pitch  plane  between  the  missile  longitudinal  (x  axis)  and  the  z 
component  of  the  total  missile  body  velocity  vector.  The  sideslip  angle 
(/))  is  the  angle  formed  in  the  yaw  plane  between  the  missile 
longitudinal  axis  and  the  y  component  of  the  total  missile  body  velocity 
vector. 

Figure  2  shows  the  inertial  reference  frame .  Note ,  that  the  2 ^ 
axis  as  shown  is  positive  down.  Shown  in  Figure  2  are  azimuth  and 
elevation  angles,  and  the  relative  range  (R)  from  the  missile  to  the 
target.  This  figure  is  primarily  used  to  describe  the  missile- target 
intercept.  The  inertial  reference  frame  can  be  transformed  into  the 
missile  body  reference  frame  by  using  Euler  angles.  The  Euler  angles  <f> , 
8 ,  and  rf>  are  used  to  describe  the  angles  of  rotation  needed  for  the 
inertial  reference  frame  to  come  into  coincidence  with  the  missile  body 
reference  frame.  There  is  a  specific  order  to  the  rotation.  First,  the 
inertial  reference  frame  is  rotated  about  the  Z axis  through  the  angle 
rj>.  Second,  the  inertial  reference  frame  is  rotated  about  the  new  Y 

i 

axis  through  the  angle  8.  Finally,  the  inertial  reference  frame  is 

rotated  about  the  new  X  axis. 

i 

3.  THE  BANK-TO-TURN  CONTROL  PROBLEM 

In  Reference  1,  Arrow  summarizes  preferred  orientation  control  for 
tactical  missiles.  Missile  steering  falls  into  two  basic  categories: 
skid-to-tum  (STT)  and  bank-to-turn  (BTT)  control.  STT  control  is 
prevalent  in  most  existing  guided  weapons.  In  STT  control  the  maneuvers 
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are  done  in  the  pitch  and  yaw  plane  with  the  missile  roll  stabilized. 

In  BTT  control  the  maneuvers  are  primarily  done  in  the  roll  and  pitch 
planes  while  control  may  be  used  to  minimize  sideslip.  BTT  control  of 
the  missile  is  similar  in  nature  to  how  pilots  might  fly  their  aircraft, 
Reference  2.  BTT  missiles  have  only  recently  begun  to  receive  more 
attention  for  their  possible  use  in  air-to-air  applications.  BTT 
vehicles  have  two  chief  advantages  over  STT  vehicles.  The  first  is  that 
BTT  vehicles  can  travel  farther  with  the  same  amount  of  fuel.  BTT 
vehicles  also  have  a  much  larger  acceleration  capability  in  the 
preferred  maneuver  plane  than  STT  vehicles. 

BTT  steering  may  be  placed  into  categories  such  as  BTT  180,  BTT  90, 
and  coordinated  BTT  (CBTT) .  BTT  180  means  the  airframe  may  be  commanded 
to  roll  as  much  as  180  degrees  before  acceleration.  BTT  90  airframes 
may  roll  up  to  90  degrees.  CBTT  implies  the  entire  maneuver  is 
coordinated  so  that  the  airframe  is  commanded  to  accelerate  in  a 
specified  direction  before  the  roll  maneuver  completes .  The  airframe 
used  in  this  study  is  a  BTT  180  airframe  and  uses  CBTT  command  logic  to 
command  the  autopilot. 

BTT  airframes  are  capable  of  high  roll  rates  which  induce 
undesirable  yawing  and  pitching  moments  in  the  airframe.  The 
asymmetrical  properties  of  BTT  airframes  lead  to  a  requirement  to 
minimize  sideslip  since  large  sideslip  angles  have  a  negative  affect  on 
pitch  and  roll  dynamics.  Another  requirement  to  minimize  sideslip  is  to 
prevent  the  air  intake  for  air  breathing  missiles  from  being  reduced  to 
the  point  that  the  propulsion  unit  fails  (this  is  called  flameout) . 

4.  MISSILE  CROSS  COUPLING 

BTT  missiles  are  asymmetrical  airframes,  and  the  lack  of  symmetry 
leads  to  increased  gyroscopic  and  inertial  coupling  of  roll,  pitch,  and 
yaw  rates.  The  roll  rate  coupling  may  be  the  most  severe  of  the  three. 
Since  BTT  control  is  used,  the  airframe  must  be  capable  of  attaining 
high  roll  rates.  These  high  roll  rates  are  coupled  into  the  pitch  and 
yaw  dynamics.  When  the  assumption  is  made  that  the  airframe  is 
symmetrical  some  of  the  cross  coupling  terms  are  assumed  to  be 
negligible  in  the  moment  equations.  This  assumption  is  primarily 
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valid  for  STT  airframes  which  usually  are  nearly  symmetrical.  This 
assumption  also  leads  to  a  simplified  nonlinear  model.  In  the  approach 
used  for  this  problem,  symmetry  is  not  assumed;  hence,  the 
cross  coupling  terms  are  not  immediately  assumed  to  be  negligible. 

Aerodynamic  cross  coupling  also  affects  the  roll,  pitch  and  yaw 
dynamics.  The  aerodynamic  flying  qualities  and  coupling  effects  are 
normally  obtained  through  wind  tunnel  testing.  The  effects  of 
aerodynamic  coupling  are  then  clearly  a  function  of  each  of  these 
coefficients  and  will  vary  from  one  airframe  to  the  next. 

5 .  AUTOPILOT 

Figure  3  shows  a  simplified  missile  system  block  diagram.  As  shown 
in  Figure  3,  the  autopilot  consists  of  the  BTT  logic  and  a  controller 
with  a  gain  scheduler.  Since  the  controller  is  based  on  a  linearized 
model,  gain  scheduling  is  required  to  implement  a  piecewise  linear 
controller.  Cain  scheduling  consists  of  gain  lookup  tables  that  allow 
gains  to  change  with  varying  parameters.  The  autopilot  takes  the 
guidance  law  commands,  which  are  commanded  accelerations  in  the  body 
frame,  and  develops  commands  for  the  actuators.  The  actuators  move  the 
missile  control  surfaces  (fins).  The  aileron,  elevator,  and  rudder 
deflection  commands  are  mixed  to  develop  commands  for  the  four 
actuators . 
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SECTION  II 
PREVIOUS  WORK 


Many  studies  for  BTT  autopilots  have  been  conducted  during  the  past 
several  years.  Approaches  to  autopilot  design  range  from  classical 
designs  to  linear  quadratic  regulator  (LQR) ,  linear  quadratic 
gaussian/loop  transfer  recovery  (LQG/LTR) ,  pole  placement,  and 
eigenstructure  assignment  designs.  An  assessment  of  BTT  autopilot 
technology  is  provided  by  Cloutier,  Evers,  and  Feely  in  Reference  3. 

1.  CLASSICAL  DESIGNS 

In  one  approach  roll,  pitch,  and  yaw  autopilots  were  independently 
designed  and  the  roll  command  from  the  roll  autopilot  fed  to  the  pitch 
and  yaw  autopilots.  Reference  4.  The  intent  was  to  coordinate  the  roll 
and  yaw  motion  and  minimize  coupling  of  roll  and  yaw  dynamics  into  the 
pitch  dynamics.  In  this  approach  the  yaw  command  follows  the  roll 
command;  this  minimizes  sideslip  and  is  termed  as  a  coordinated 
bank-to-tum  steering. 

In  another  design,  a  dither  adaptive  approach  was  used  to 
independently  design  roll,  pitch  and  yaw  autopilots.  Reference  5.  In 
this  approach  a  dither  signal  is  injected  and  the  adaptive  network  is 
used  to  tune  the  rate  and  acceleration  loop  gains  of  the  autopilots . 

2.  MODERN  DESIGNS 

An  LQG/LTR  pitch/yaw  autopilot  based  on  a  model  with  the  roll 
dynamics  decoupled  from  the  pitch/yaw  dynamics  was  designed  with  pole 
placement  techniques  used  for  the  roll  autopilot,  Reference  6. 

Gyroscopic  coupling  of  the  roll  dynamics  into  the  pitch  dynamics  was 
included  by  developing  linear  pitch/yaw  models  at  fixed,  nonzero  roll 
rates.  The  autopilot  gains  were  scheduled  as  a  function  of  roll  rate 
and  dynamic  pressure. 

In  another  variation,  an  LQR  design  with  the  pitch  dynamics 
decoupled  from  the  roll/yaw  dynamics  was  designed,  Reference  7.  Both 
the  pitch  and  roll/yaw  autopilots  were  developed  with  LQR  techniques. 
Pitch  dynamics  entered  the  roll/yaw  design  since  the  roll/yaw  model  was 
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linearized  about  fixed  angles -of -attack.  The  design  was  digital  and 
included  a  computational  delay.  The  resulting  controller  was  scheduled 
as  a  function  of  angle -of- attack.  Dynamic  pressure  was  excluded  from 
the  schedule  since  performance  comparisons  were  done  at  fixed  altitudes. 
Much  of  this  work  is  used  as  a  basis  for  this  thesis. 

In  a  comparison  of  a  modern  control  approach  to  a  classical 
approach  Arrow  and  Williams  state  that  classical  approaches  can  fall 
short  of  modern  approaches  in  ease  of  finding  a  control  law  for  a  highly 
coupled  system,  Reference  8.  They  also  state  the  classical  approach 
deals  better  with  unmodeled  constraints. 

In  Reference  9,  Nesline  and  Zarchan  demonstrate  the  benefits  of 
analyzing  modern  control  designs  using  classical  frequency  response 
techniques.  The  purpose  is  to  identify  unmodeled  disturbances  that 
potentially  may  destabilize  the  modern  design. 

3.  DECOUPLED  MISSILE  DYNAMICS 

The  thrust  of  these  studies  still  centers  around  decoupled  designs. 
In  these  decoupled  approaches  the  nonlinear  model  is  simplified  before 
the  linear  model  is  developed.  In  most  approaches  to  simplifying  the 
nonlinear  model  the  basic  assumption  is  that  the  effects  of  dropping  the 
inertial  cross  coupling  terms  are  minimal.  In  most  cases  this  may  in 
fact  be  true  since  aerodynamic  forces  and  moments  tend  to  dominate  the 
missile  dynamics.  The  validity  of  the  assumption  is  tested  by 
evaluating  the  performance  of  the  autopilot  in  a  detailed  nonlinear 
missile  simulation  which  includes  all  dynamics. 

Multiple  input,  multiple  output  designs  are  available  with  modern 
control  approaches,  which  leads  one  to  expect  autopilot  designs  based  on 
the  full  roll/pitch/yaw  missile  dynamics  (roll/pitch/yaw  autopilots) 
should  be  more  prevalent  when  modern  control  approaches  are  used.  One 
reason  for  the  lack  of  roll/pitch/yaw  designs  may  be  the  influence  of 
past  classical  autopilot  design  approaches.  Previously,  most  air-to-air 
missiles  were  STT  and  with  symmetrical  qualities.  This  simplified  the 
design  considerations,  and  allowed  treatment  of  the  problem  with 
classical  designs  for  single  input  single  output  systems.  Also, 
depending  on  the  airframe,  a  decoupled  autopilot  design  may  be  desirable 
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since  the  decision  to  decouple  the  design  is  based  on  some  insight  into 
the  airframe  characteristics.  Based  on  the  linear  model  for  the 
airframe  used  in  this  study  it  still  appears  that  a  roll/pitch/yaw 
design  based  on  the  full  nonlinear  model  should  lead  to  increased 
robustness  and  at  the  very  least  offer  a  means  of  evaluating  the 
decoupled  approaches . 
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SECTION  III 


NONLINEAR  AIRFRAME 


1.  DESCRIPTION 

The  extended  medium  range  air-to-air  technology  (EMRAAT)  airframe 
is  used  in  this  study.  The  EMRAAT  is  a  BTT  airframe  with  four  control 
surfaces  (fins)  that  are  controlled  by  actuators.  Figure  4  shows  a 
diagram  of  the  airframe  used.  As  can  be  seen  in  this  figure  the 
airframe  is  asymmetric.  The  top  of  the  airframe  is  flat  and  the  tail 
fins  are  not  cruciform.  Cruciform  fins  form  four  90  degree  angles 
unlike  the  fins  shown  in  Figure  4.  This  airframe  has  one  preferred 
maneuver  plane  which  is  normal  to  the  bottom  surface  of  the  missile. 
This  implies  the  airframe  may  have  to  roll  as  much  as  180  degrees  to 
achieve  a  particular  maneuver. 

2.  EULER  EQUATIONS  OF  MOTION 

Missiles,  as  with  any  airframes,  have  six  degrees  of  freedom  which 
are  described  by  the  Euler  equations  of  aerodynamic  motion.  The  Euler 
equations  describe  the  forces  and  moments  acting  on  the  missile.  The 
six  degrees  of  freedom  are  based  on  a  missile  body  coordinate  frame  as 
shown  in  Figure  5 .  Three  equations  describe  the  forces  acting  on  the 
missile  body.  The  three  force  equations  are  as  follows: 


F  +  T  +  W 


F  +  T  +  W 
y  y  y 


F  +  T  +  W 
z  z  z 


^u  +  qw  -  rvj 
~^v  +  ru  -  pwj 


jj[w  -  qu  +  pvj 


(1) 


(2) 


(3) 
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Top  View 
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Side  View 


Figure  4.  EMRAAT  Airframe 
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In  Equations  (1),  (2),  and  (3),  u,  v,  and  w  denote  the  missile  body 

velocities.  The  variables  p,  q,  and  r  denote  the  respective  roll,  pitch 

and  yaw  rates.  F  ,  F  ,  and  F  are  force  components  that  are  defined  by 
x  y  z 

data  for  each  specific  airframe.  Wind  tunnel  testing  is  used  to  obtain 

this  data.  T  ,  T  ,  and  T  are  forces  as  a  result  of  the  thrust 
x  y  z 

components  acting  on  the  missile.  This  information  is  based  on  the 
motor  used  for  the  airframe.  W  ,  W  ,  and  W  are  the  forces  acting  on 
the  missile  due  to  the  missile  weight  components.  W  is  the  missile 
weight  and  g  is  gravity. 

Three  equations  describe  the  moments  acting  on  the  missile.  These 
three  moment  equations  are  as  follows 


1  - 

I  p  - 
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f1  -  I  1 
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^  zz  xxj 
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I  r  - 
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f1  -  1  1 

t  xx  yyJ 

pq  +  i„(q2 

■A 

•‘..l 
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•  p) 
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(6) 


Where  the  moments  of  inertia  are 
products  of  ineLtia  are  given  by 
and  (6)  can  be  written  in  matrix 


given  by  I 


xy  xz’ 
form 


xx  n 

and  I  . 
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and  I  The 
zz 

Equations  (4),  (5) 
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Where 
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(I  -  I  )qr  +  I  (r2  -  q2)  -  I  pq  +  I  rp 
yy  zz  ^  yz  n  xzrn  xy  r 


(I  -  I  )rp  +  I  (p2  -  r2)  -  I  qr  +  I  pq 
zz  xx  r  XZ  r  XYn  yz  ^ 


2  2 

(I  -  I  )pq  +  I  (q  -  p)-I  rp  +  I  qr 

XX  YY  XY  ^  r  YZ  r  XZn 


(9) 


3.  NONLINEAR  AIRFRAME  STATE  EQUATIONS 

This  section  presents  the  development  of  the  nonlinear  airframe 
model  as  used  in  this  study  and  is  based  on  the  aerodynamic  data  for  the 
airframe.  This  development  leads  to  the  five  nonlinear  state  equations 
which  are  later  linearized  to  get  the  linear  model  for  the  autopilot 
design.  The  force  equations  are  developed  first  and  then  the  moment 
equations  are  developed. 

As  was  done  in  Reference  7,  and  other  autopilot  designs,  Equation 
(1)  is  not  included  as  part  of  the  nonlinear  airframe  state  equations. 
Equation  (1)  is  included  in  the  nonlinear  six-degree-of-freedom  missile 
simulation.  For  this  study  thrust  is  assumed  to  be  zero,  so 


T  -  T  -  T  -  0 
x  y  z 


(10) 
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The  force  components  for  weight  are  given  as 


Wy  -  Wcos(0 )sin(^) 

(ID 

Wz  -  Wcos(0 )cos(^) 

(12) 

Where  9  and  ^  are  two  of  the  three  Euler  angles  that  describe  the 
inertial  orientation  of  the  missile  body  with  respect  to  the  body 
reference  frame  (Figure  5).  Substituting  Equations  (11)  and  (12)  into 
Equations  (2)  and  (3)  and  solving  for  F  and  F  yields 

y  z 


F^  -  (v  +  ru  -  pw)  -  Wcos(0 )sin(^) 


W  • 

F  -  —  (w  -  qu  +  pv)  -  W cos(9)cos(<j>) 
z  g 


The  angle -of  attack  and  sideslip  angle  are  given  by 


a  -  tan 


l(5) 


P  —  tan 


im 


(13) 

(14) 


(15) 


(16) 


During  flight  a  and  P  will  be  small,  so  the  following  small  angle 
approximations  can  be  used: 


a 


(17) 


(18) 
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The  total  missile  velocity  (V)  is  given  by 


V 


1/2 

[(U2  +  v2  +  w2)j 


(19) 


Since  the  missile  speed  is  nearly  constant  and  the  largest  component  of 
the  velocity  is  along  the  missile's  x  axis  the  following  approximations 
can  be  used 

u  =  V  (20) 

u  =  0  (21) 

It  follows  from  Equations  (17)  and  (18)  that 


(22) 

(23) 


Using  these  approximations,  Equations  (13)  and  (14)  can  be  rewritten  to 
eliminate  u,  v,  and  w,  and  introduce  a,  0,  and  V,  as  follows: 


F  -  —  03  +  r  -  pa)  -  Wcos(0 )sin(*)  (24) 

y  s 


F  -  (a  -  q  +  v0)  -  Wcos(0)cos(^)  (25) 

z  g 
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The  Appendix  lists  the  aerodynamic  data  for  the  EMRAAT  airframe  for 
two  flight  conditions,  the  first  condition  is  40,000  feet  at  mach  2.5, 
and  the  second  condition  is  at  65,000  feet,  at  mach  3.0.  The  Appendix 
provides  the  aerodynamic  force  derivatives  which  are  used  to  specify  the 
left  hand  sides  of  Equations  (24)  and  (25).  Equations  (24)  and  (25)  can 
be  written  as  (Reference  10) 


Fy  -  (QS) 
Fz  -  (QS) 


\C  f3  +  C  p  +  C  r  +  C  6  +  C  S 

!  Y  “  Y  YpYr 

^  $  p  t  S  p  4  r 


[c  a  +  C  a  +  C  q  +  C  S  I 

8  S.  N  ^  N  q  I 

v  a  a  q  4q 


) 


(26) 


(27) 


Substituting  Equation  (26)  and  (27)  into  Equations  (24)  and  (25)  and 
solving  for  a  and  p  leads  to 


-  1 


a  — 
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1  .  - - 
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q  s  q  J  v 


(28) 
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fi 


gQS  f 

pot  -  r  +  -  C  fl  +  C  p+C  r  +  C  S 
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WV  v  0  P  I  ip 


g 

+  -  cos(6)cos(<fi) 
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The  state  Equations  (28)  and  (29)  are  the  two  force  equations  used 
in  development  of  the  linear  model.  The  next  step  is  the  development  of 
the  state  equations  for  the  moments.  Equation  (7)  can  be  solved  to 
yield 
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-  j 
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-  j'Vp.q.r) 


(30) 


Where  H(p,q,r)  is  given  by  Equation  (9).  To  complete  Equation 
(30),  the  aerodynamic  moments  1,  m,  and  n  are  provided  as  part  of  the 
data  from  the  aerodynamic  tables  supplied  in  the  appendix  and  can  be 
written  in  the  form  (Reference  10) 
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m 
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(32) 


(33) 
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The  nonlinear  model  development  is  now  complete.  The  five  state 
nonlinear  model  from  which  the  linearized  model  is  developed  is 
described  by  Equations  (28) ,  (29) ,  and  (30) .  The  Appendix  lists  the 
values  for  the  variables  used  in  the  equations.  The  states  used  in  the 

model  are 


x 


[■ 


P  p  q 


(34) 


4.  MODEL  SIMPLIFICATIONS  NOT  USED 

A  closer  look  at  Equations  (28)  through  (33)  reveals  the 
aerodynamic  and  inertial  cross  coupling  present  in  a  nonlinear  airframe. 
The  model  as  it  was  presented  is  the  model  that  will  be  used  for  the 
autopilot  design.  It  should  be  noted  that  in  other  approaches  the  model 
dynamics  were  simplified  which  leads  to  reduced  order  linearized  models. 
Some  assumptions  that  can  be  used  to  reduce  the  order  of  the  model  so 
the  dynamics  can  be  separated  into  a  pitch  model  and  a  roll/yaw  model 
are  presented  here. 

In  terms  of  the  aerodynamic  data  provided  for  the  EMRAAT  airframe, 
Equations  (26),  (27),  (31),  (32),  and  (33)  have  no  terms  that  couple  the 
pitch  dynamics  to  the  yaw  or  roll  dynamics.  These  equations  do  Indicate 
there  is  aerodynamic  coupling  between  the  yaw  and  roll  dynamics.  This 
partially  supports  using  two  simplified  models,  one  for  the  pitch 
dynamics  and  another  for  the  roll/yaw  dynamics. 

The  force  equations  are  simplified  first.  Equation  (28)  shows  some 
cross  coupling  of  roll  and  yaw  dynamics  into  the  pitch  dynamics. 

Equation  (29)  shows  cross  coupling  of  pitch  and  roll  dynamics  into  the 
yaw  dynamics.  To  remove  the  roll/yaw  dynamics  from  Equation  (28)  the 
term  p£  can  be  assumed  to  be  negligible  since  0  is  assumed  to  be  near 
zero.  To  remove  the  pitch  dynamics  from  Equation  (29)  the  term  pa  can 
be  changed  to  pa^  where  is  a  constant  based  on  the  assumption  that  a 
is  slowly  changing. 
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The  moment  equations  are  now  simplified.  Based  on  the  assumption 
that  the  airframe  is  nearly  symmetrical,  the  moments  of  inertia  I  and 
I  are  treated  as  being  equal  and  the  products  of  inertia  1,1,  and 

ZZ  XY  XZ 

I  are  neglected.  This  assumption  reduces  the  inertia  matrix  J  in 
Equation  (8)  to 


I 

xx 


J  - 


0 


0 


0 


(35) 


Equation  (9)  may  also  be  reduced  to 
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(I  -I  )rp 
zz  xx7  r 


-(I  -I  )pq 

XX 


(36) 


Equation  (30)  when  used  with  Equations  (35)  and  (36)  can  now  be 
separated  into  a  model  for  the  pitch  dynamics  and  a  model  for  the 
roll/yaw  dynamics.  This  is  one  possible  approach  to  simplifying  the 
model.  Other  approaches  have  separated  the  model  into  roll  dynamics  and 
pitch/yaw  dynamics.  Classical  approaches  have  separated  the  model  into 
roll,  pitch,  and  yaw  models  and  later  added  feedback  paths  between 
controllers  to  account  for  cross -coupling. 
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5.  MISSILE  CONTROL  SURFACES 

Control  for  the  missile  accelerations  along  the  y  and  z  axis  and 
the  missile  body  accelerations  is  provided  by  the  missile  control 
surfaces .  The  control  surfaces  are  shown  in  Figure  6 .  As  can  be  seen 
in  Figure  6  the  tail  configuration  is  not  symmetrical.  Clearly, 
controls  will  also  interact  to  compound  the  overall  problem  of  coupling. 
Equation  (37)  is  the  mixing  logic  required  to  transform  the  four  control 
surfaces  shown  in  Figure  6  to  the  aileron,  elevator  and  rudder  controls 
(5,5,  and  5  ).  Although  not  actual  control  surfaces,  5,5,  and  5 

p  q  r  p  q  t 

are  used  in  the  aerodynamic  data  provided  the  EMRAAT  airframe.  Another 

reason  5,5,  and  5  are  used  is  that  both  aerodynamic is ts  and  flight 
p  q  r 

controls  engineers  find  it  easier  to  think  of  the  control  surfaces  in 
terms  of  aileron,  elevator,  and  rudder. 
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6.  SIX-DEGREE-OF-FREEDOM  (6-DOF)  MISSILE  SIMULATION 

The  6-DOF  simulation  embodies  the  full  nonlinear  dynamics  for  a 
guided  missile  system.  Figure  7  shows  a  simplified  subsystems  model  for 
a  guided  missile.  These  subsystems  are  implemented  in  the  6-DOF 
simulation  with  considerably  more  detail  than  can  be  shown  on  this 
diagram.  The  bold  lines  on  the  diagram  indicate  signal  flow.  The 
dashed  lines  indicate  measurements. 

The  missile  guidance  consists  of  the  seeker  and  a  guidance  law. 

The  seeker  modeled  is  an  active  seeker  which  means  it  is  emitting  energy 
to  scan  the  target.  Gimbal  dynamics  are  not  included  in  the  seeker 
model.  The  information  provided  by  the  seeker  to  the  guidance  law 
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consists  of  the  missile- target  line-of-sight  (LOS)  angles,  LOS  angle 
rates,  range,  and  range  rate.  The  guidance  law  is  proportional 
navigation  and  provides  lateral  acceleration  commands  to  the  BTT  logic. 

The  inertial  measurement  unit  (IMU)  consists  of  accelerometers  and 
rate-gyros  to  measure  the  missile  body  accelerations  and  missile  body 
roll,  pitch,  and  yaw  rates  (turning  rates).  Prefilters  are  used  to 
filter  the  measurements.  The  parameters  are  used  throughout  the 
simulation. 

The  BTT  logic  uses  the  lateral  acceleration  commands  from  the 
guidance  law  and  the  IMU  measurements  to  form  the  autopilot  commands. 
Included  in  the  BTT  logic  is  an  estimator  to  estimate  a  and  /?  in  the 
presence  of  wind  and  sensor  noise. 

The  autopilot  is  implemented  in  the  simulation  as  a  discrete 
controller  with  a  computational  delay  included.  The  autopilot  gains  are 
scheduled  to  develop  a  piece  wise  linear  autopilot. 

The  actuators  receive  their  commands  from  the  autopilot  to  control 
the  total  airframe  response.  In  the  simulation,  each  actuator  is 
modeled  as  a  fifth  order  electromechanical  model  and  includes  rate 
limiting.  The  actuators  move  the  missile  control  surfaces.  The 
autopilot  commands  to  the  actuator  are  delayed  by  6  ms  to  simulate  a 
computational  delay. 

In  Figure  7,  the  block  containing  the  airframe  dynamics  includes 
the  full  nonlinear  missile  dynamics  for  the  EMRAAT  airframe.  This  block 
generates  the  airframe  response  to  the  actuator  inputs. 

Also  included  in  the  simulation  is  a  maneuvering  target  capable  of 
performing  9g  maneuvers  to  evade  the  missile.  The  simulation  calculates 
the  time  left  for  intercept  to  occur  and  if  this  time  reaches  a 
specified  minimum  threshold  the  target  performs  an  evasive  maneuver. 

Two  feedback  loops  should  be  apparent  in  Figure  7.  The  inner  loop 
consists  of  the  autopilot,  BTT  logic,  actuators,  airframe  dynamics  and 
the  IMU.  The  outer  loop  includes  the  guidance  and  the  inner  loop. 
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Figure  6.  Control  Surfaces 
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AUTOPILOT 


Figure  7.  Guided  Missile  Subsystems 


SECTION  IV 

LINEAR  MODEL  DEVELOPt-ENT 

1.  DEVELOPMENT  OF  THE  LINEAR  MODEL 

The  five  state  nonlinear  model  of  Section  III  is  now  used  to 
develop  the  linearized  model.  The  nonlinear  system 

x(t)  -  f (x(t) ,u(t))  (38) 

is  linearized  by  using  Reference  11  and  leads  to  the  linearized  model 

x(t)  -  Ax(t)  +  Bu(t)  (39) 


Where  A  is  given  by 


A 


df 

dx  ^ 

nominal 


(40) 


and  B  is  given  by 


B 


- 

du'x 

nominal 


(41) 


and  where  Ax(t)  and  Au(t)  have  been  redefined  to  be  x(t)  and  u(t)  The 
system  state  and  control  vectors  are 

T 

x(t)  -  |a  p  p  q  rj  (42) 

and 

T 

u(t)  -  6^  S  J  (43) 
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Using  Equations  (40)  and  (41) ,  the  model  is  now  ready  to  be  linearized 
to  yield  the  necessary  equations  for  each  element  of  the  linear  A  and  B 
system  matrices.  To  simplify  the  equations  let 


and 
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*SC.. 
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wv 


(44) 
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Equations  (28)  and  (29)  are  first  linearized  to  yield  the  first  two  rows 
of  A  and  B.  The  first  element  of  A  is  found  by  taking  the  partial 
derivative  of  (28)  with  respect  to  a  and  leads  to 
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(46) 


Taking  the  partial  derivative  of  Equation  (28)  with  respect  to  0  leads 
to 
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12 


Z 
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(47) 


Taking  the  partial  derivative  of  Equation  (28)  with  respect  to  p  leads 
to 
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Taking  the  partial  derivative  of  Equation  (28)  with  respect  to  q  leads 
to 


27 


a 


14 


(1  +  CK  Q~  ) 
q 


Z 

d 


(49) 


Taking  the  partial  derivative  of  Equation  (28)  with  respect  to  r  leads 
to 


15 


-  0 


(50) 


The  first  element  in  B  is  found  by  taking  the  partial  derivative  of  (28) 
with  respect  to  5  and  leads  to 
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The  partial  derivative  of  Equation  (28)  with  respect  to  S  leads  to 
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and  the  partial  derivative  of  Equation  (28)  with  respect  to  S  leads  to 


b  -  0  (53) 
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To  find  the  second  row  elements  for  A  and  B,  the  procedure  to  find  the 
row  one  elements  is  repeated  with  the  exception  that  Equation  (28)  is 
now  replaced  with  Equation  (29)  and  the  partial  derivatives  are  taken 
with  respect  to  q.  The  row  2  elements  for  A  and  B  are  found  to  be 
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Extra  steps  are  necessary  to  linearize  the  moment  equations  which 
will  lead  to  rows  3,  4  and  5  for  A  and  B.  The  first  step  is  to  multiply 
Equation  (30)  by  the  inertia  matrix  J  and  then  linearize  the  right  hand 
side  of  the  result.  This  can  be  done  since  J  is  independent  of  the 
states.  J  will  be  included  later  in  the  linear  model  development.  This 
leads  to 
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Linearization  of  Equation  (62)  leads  to  the  matrices  M  and  N  which  will 
become  part  of  matrices  A  and  B.  The  matrix  M  is  given  by 
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For  m  to  m  .  and  n  to  n  the  a  derivatives  must  be  included. 
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Taking  the  partial  derivatives  of  row  2  of  Equation  (62)  with  respect  to 
each  of  the  elements  in  x  and  u  leads  to 
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Finally,  taking  the  partial  derivatives  of  row  3  of  Equation  (62)  with 
respect  to  each  of  the  elements  in  x  and  u  leads  to 
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Since  the  inertia  matrix  (J)  is  independent  of  the  states  and 
controls,  the  matrices  M  and  N  can  now  be  multiplied  by  J  \  The  system 
matrices  A  and  B  can  be  written  as 
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Equations  (89)  and  (90)  are  in  the  form  used  to  compute  the 
linearized  models  for  the  selected  design  points. 

2.  LINEAR  MODEL  FOR  SYSTEM  1 

The  linearized  model  for  an  altitude  of  40,000  feet,  mach  2. 
angle -of- attack  of  10  degrees  is  given  by 
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-150.58 

(92) 


Sideslip  angle,  and  body  rates  are  all  zero  for  this  model.  For 
purposes  of  example  this  model  will  be  called  linear  system  number  1. 

3.  LINEAR  MODEL  FOR  SYSTEM  2 

The  linearized  model  for  an  altitude  of  65,000  feet,  mach  3.0,  and 
angle - of- attack  of  10  degrees  is  given  by 
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Sideslip  angle,  and  body  rates  are  all  zero  for  this  model.  For 
purposes  of  example  this  model  will  be  called  linear  system  number  2. 
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SECTION  V 

DESIGN  CONSIDERATIONS 


A  number  of  consideraCions  play  an  Important  role  during  autopilot 
design.  Angle  limitations  are  placed  upon  the  sideslip  angle  and 
angle -of -attack.  The  design  must  take  into  account  the  method  for 
commanding  the  airframe  such  as  BTT-180  or  BTT-90.  The  actuators  must 
also  be  considered  since  they  too  impact  autopilot  performance. 


1.  SIDESLIP  ANGLE 

The  sideslip  angle,  f),  must  be  limited  to  ±5  degrees  to  limit  the 
coupling  into  the  roll  and  pitch  channels.  The  commanded  sideslip  from 
the  BTT  logic  is  zero.  The  EMRAAT  airframe  has  a  strong  tendency  to 
return  to  zero  (or  near  zero)  sideslip  so  this  condition  is  easily  met. 
Commands  for  small  sideslip  angles  may  be  desirable  in  terms  of  reducing 
undesirable  effects  from  seeker-roll  channel  coupling,  but  this  effect 
is  not  addressed  in  this  project. 

2.  ANGLE-OF-ATTACK 

The  angle -of -attack,  a,  must  be  limited  from  between  -5  degrees  to 
+25  degrees.  This  is  the  range  for  which  the  aerodynamic  data  provided 
for  the  airframe  is  valid.  Because  of  the  airframe  design  and  control 
surface  configuration,  negative  angles -of -attack  reduce  the  control 
effectiveness.  The  BTT  logic  therefore  commands  angles -of -attack  from 
between  0  to  +20  degrees. 

3 .  BTT  LOGIC 

The  BTT  logic  determines  the  manner  in  which  the  airframe  should 
maneuver  in  response  to  guidance  law  commands.  The  BTT  logic  used  in 
this  study  is  designed  to  command  the  missile  to  roll  about  the  velocity 
vector.  To  do  this  rate  commands  for  roll,  pitch,  and  yaw  rates  and 
angle  commands  for  angle -of- attack  and  sideslip  are  provided  to  the 
autopilot.  Using  this  BTT  logic  the  response  times  for  both  the  roll 
rate  and  angle  of  attack  become  important. 
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4.  ACTUATOR  MODEL 

The  actuator  model,  for  each  of  the  control  surfaces,  is  a  fifth 
order  electro -mechanical  model.  Since  three  actuators  are  used,  they 
comprise  a  15  state  actuator  system.  This  system  is  left  out  of  the 
autopilot  design.  It  is  important  that  the  bandwidth  of  the  actuators 
be  as  high  as  possible,  to  accommodate  the  autopilot,  while  not  being  so 
high  that  rate  limiting  occurs  within  the  actuators.  The  bandwidth  for 
these  actuators  is  27  Hz. 

5.  ADDITIONAL  CONSIDERATIONS 

This  study  addresses  only  the  missile  rigid  body.  If  the  problem 
of  the  missile  flexible  modes  were  to  be  addressed  then  bandwidth 
restrictions  on  the  controller  become  more  severe  since  it  is  desirable 
to  attenuate  these  modes . 

The  design  is  done  for  the  case  when  the  fuel  has  been  exhausted. 
This  simplifies  the  approach  in  that  changing  missile  mass  and  the 
changing  center  of  mass  are  not  considered.  All  simulations  of  the 
design  assume  near- intercept  conditions  which  also  means  the  fuel  supply 
is  exhausted  (or  nearly  exhausted) .  This  is  a  good  point  to  start  the 
design  since  it  is  in  near- intercept  conditions  where  the  autopilot 
performance  is  most  critical.  One  of  the  next  steps  in  the  design  is  to 
include  the  case  of  a  thrusting  motor  with  changing  missile  mass.  The 
thrusting  motor  is  not  within  the  scope  of  this  thesis. 

The  design  must  address  the  fact  that  the  autopilot  is  a  digital 
implementation  in  a  microprocessor.  The  autopilot  receives  discrete 
commands  from  the  guidance  law  at  a  rate  of  100  Hz  and  provides  discrete 
commands  for  the  actuators  at  the  same  rate.  Also  of  concern  is  the 
computational  delay  which  is  from  the  time  the  guidance  law  commands  the 
autopilot  to  the  time  the  autopilot  commands  the  actuators.  The 
estimated  computational  delay  is  6  mS. 
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SECTION  VI 


APPROACH 


1.  LINEAR  REGULATOR  MODEL 

Figure  8  shows  a  simplified  model  of  the  missile  control  system. 

It  should  be  noted  that  the  guidance  law,  BTT  logic,  and  seeker  require 
more  information  about  the  missile  dynamics  than  is  shown  in  Figure  8. 
The  figure  is  intended  to  show  the  interrelationship  of  the  controller 
to  the  other  systems  and  the  airframe  dynamics.  The  autopilot  consists 
of  the  BTT  logic  and  controller  (with  a  gain  schedule  implemented) . 

The  function  of  the  BTT  logic  is  to  transform  the  lateral 
acceleration  commands  from  the  guidance  law  to  commands  that  relate  to 
controlling  the  airframe.  In  this  approach  the  BTT  logic  generates 
angle  rate  commands  for  roll,  pitch  and  yaw,  and  angle  commands  for 
angle -of- attack  and  sideslip.  The  BTT  logic  diagram  is  shown  in  Figure 
9  (Reference  7).  The  BTT  logic  works  well  for  this  problem,  so  no  work 
is  done  with  this. 

The  controller  is  designed  using  linearized  quadratic  regulator 
(LQR)  methods  and  is  based  on  the  linearized  model  from  Section  IV  which 
is  linearized  about  specific  flight  conditions  and  then  discretized 
about  a  sample  rate  of  10  ms.  The  design  also  takes  into  account  a 
computational  delay  of  6  ms.  A  simplified  diagram  of  the  regulator  is 
shown  in  Figure  10.  The  rest  of  this  section  develops  the  regulator 
shown  in  Figure  10.  Since  the  true  system  model  is  nonlinear  and  the 
flight  conditions  vary  constantly,  gain  scheduling  is  implemented  to 
develop  a  piecewise  linear  regulator. 
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Figure  8.  Missile  System  Block  Diagram 
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Figure  9.  BTT  Logic 


2.  LINEAR  QUADRATIC  REGULATOR  DESIGN 

Of  the  five  states  in  x,  a,  fi  and  p  require  more  precise  control 
than  q  and  r.  As  was  previously  explained  in  Section  V,  there  is  no 
problem  in  controlling  /9.  The  linear  model  is  augmented  with  integral 
states  for  a  and  p  to  minimize  their  steady- state  errors  and  improve 
command  tracking.  The  augmented  system  can  be  specified  as: 

T 

xt  -  [«  |  Pj]  (95) 

Where  and  p^  are  a  time  integration  of  the  states  a  and  p.  The 
System  matrices  A  and  B  are  augmented  with  zeros  and  ones  to  include  the 
new  states.  This  leads  to 


(96) 
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In  the  LQR  approach  the  optimal  control  for  u  is: 

u  -  -  K  x  (98) 
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Where  K  is  the  gain  matrix.  Transmission  zero  weighting  is  used  to 
specify  the  quadratic  weighting  matrices  Q  and  R  for  the  performance 
index  (Reference  12) : 


J 
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uTR 


■) 


dt 
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In  this  approach,  the  same  number  of  outputs  are  specified  as  are 
inputs.  Transmission  zeros  are  specified  in  a  manner  to  introduce 
performance  criteria.  In  the  resulting  closed  loop  system  as  the  loop 
gain  increases  from  a  very  small  value,  the  eigenvalues  for  the  closed 
loop  system  migrate  from  the  open  loop  system  eigenvalues  to  the 
transmission  zeros.  Following  is  a  description  of  the  approach  used  in 
Reference  7,  but  modified  for  this  problem.  The  advantage  of  this 
approach  is  the  ability  to  specify,  to  some  extent,  the  structure  of  the 
the  closed  loop  system  eigenvalues  using  classical  performance 
parameters  such  as  damping  factors  and  natural  frequencies.  In 
Reference  7,  two  regulators  were  designed,  one  for  the  pitch  autopilot 
and  one  for  the  yaw/roll  autopilot.  The  approach  is  modified  for  one 
regulator  design  for  a  roll/pitch/yaw  autopilot. 

For  the  state  vector 

T 

-  [<*  P  p  q  r  Oj  pj  (ioo) 

the  output  vector  y  is  selected  to  be 

T 

y  “  [y.  y,  yp]  (101) 

Where  the  elements  y  ,  y  ,  and  y  are  determined  in  accordance  with  the 

•  b  p 

autopilot  designs  in  Reference  7.  The  same  numbers  are  used  for  y  ,  y  , 
and  y  as  in  Reference  7  since  the  autopilots  in  Reference  7  are  used 

p 
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for  comparison  of  performance.  The  linear  system  given  in  Equations 
(91)  and  (92)  is  used  for  this  description.  The  first  element  in 
Equation  (101)  to  be  specified  is  y  .  From  Reference  7  this  first 

<X 

starts  in  the  s-plane,  and  is  written  as 
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(102) 


In  the  time  domain 


y  -a  +  aa  +  ba  (103) 

a  a  a  I 


Next  a  is  found  from  the  linearised  state  equations  to  be 

a  -  -1.6a  +  q  +  -  0.15  (104) 

<i 

Where  the  coefficients  of  Equation  (104)  are  the  coefficients  of  the 
linear  system  matrices  A  and  B  in  Equations  (91)  and  (92) .  The 
coefficient  for  the  control  5^  is  neglected  since  it  is  small.  Equation 
(104)  is  now  substituted  into  Equation  (103)  to  yield 


y  -  (-1.6  +  a  )a  +  q  +  b  a  (105) 

a  a  a  I 


In  Reference  7,  Y  (s)  is  specified  as 

b 

Y  (s)  -  (s  +  a  )/9(s) 
P  P 


(106) 


44 


In  the  time  domain 


* +  V 


(107) 


Where  from  the  linearized  state  equations 


0  -  -0.50  +  0.2p  -  r 


(108) 


It  follows  that 


y„  -  (-0.5  +  a  )0  +  0.2p  -  r  (109) 

P  P 

As  done  in  Reference  7,  the  last  element  of  Equation  (101)  is  specified 
as 
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In  the  time  domain: 


y  ■  P  +  a  P  (111) 

p  pi 

Equations  (105) ,  (109)  and  (111)  now  specify  the  elements  of  y  in 
Equation  (101).  To  find  the  weighting  matrix  Q,  the  output  is  defined 
as: 


y  -  Cx 


(112) 
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Next  define  Q  to  equal  0.01CTC  so  that 

yTy  -  0.01xTCTCx  -  xTQx  (113) 

Where  CTC  is  scaled  by  0.01  so  it  will  be  in  the  same  numerical  range  as 
the  control  weighting  matrix  R,  this  leads  to 


Q  -  0.01CTC  (114) 

The  control  weighting  matrix  R  started  with  the  control  weighting 
matrices  used  in  Reference  7.  In  Reference  7  the  non- zero  off  diagonal 
elements  in  R  were  chosen  to  minimize  control  surface  coupling.  By 
trial  and  error,  the  diagonal  elements  of  R  were  varied  from  those  used 
in  Reference  7  to  achieve  good  performance  with  the  autopilot  operating 
in  the  presence  of  noisy  measurements.  This  choice  leads  to  less  than 
optimal  deterministic  performance.  The  matrix  R  as  used  for  this 
problem  is: 


10.2 

0 

1.8 

0 

2.5 

0 

1.8 

0 

9.2 

(115) 


Table  1  lists  the  transmission  zeros  used  for  this  design  and  specifies 
the  coefficients  for  Equations  (105),  (109)  and  (111). 
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TABLE  1.  TRANSMISSION  ZEROS  USED  FOR  LINEAR  SYSTEM  1 


COEFFICIENTS (S) 

VALUE ( S ) 

TRANSMISSION  ZERO(S) 

a  and  b 

11.7  and  39.7 

s  -  -  5.9  ±  j 2 . 3 

at  at 

a 

$ 

5.5 

s  -  -  5.5 

a 

5.0 

s  -  -  5. 0 

p 

Using  the  data  from  Table  1  in  Equations  (105),  (109),  (111)  and  using 


Equation  (114) 

to  scale  Q, 

leads 

to: 

1.0 

0.0 

0.0 

0.1 

0.0 

4.0 

0.0 

0.0 

0.3 

0.0 

0.0 

-0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.1 

Q  - 

0.1 

0.0 

0.0 

0.0 

0.0 

0.4 

0.0 

0.0 

-0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

4.0 

0.0 

0.0 

0.4 

0.0 

15.7 

0.0 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

0.3 

Table  2  lists  the  eigenvalues  for  the  open  and  closed  loop  system 


based  on  an  LQR  design  for  linear  system  number  1  (the  linear  model  in 
Figure  8) .  Also  included  in  Table  2  are  eigenvalues  for  the  closed  loop 
system  for  a  large  loop  gain  to  show  that  in  fact  as  the  loop  gain  is 
increased  the  eigenvalues  do  migrate  to  the  transmission  zeros.  The 
large  loop  gain  is  accomplished  by  reducing  R  by  a  large  scale  factor. 
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TABLE  2.  OPEN  AND  CLOSED  LOOP  EIGENVALUES 
FOR  LINEAR  SYSTEM  NUMBER  1 


OPEN  LOOP 

CLOSED  LOOP 

LARGE  LOOP  GAIN 

0.0  +  j  0.0 

-  2.8  +  j  0.0 

-5.0  +  j  0.0 

0.0  +  j  0.0 

-  5.1  +  j  0.0 

-5.5  +  j  0.0 

-0.6  +  j  0.0 

-  5.0  +  j  9.0 

-5.8  +  j  2.4 

-1.1  +  j  9.1 

-  5.0  -  j  9.0 

-5.8  -  j  2.4 

-1.1  -  j  9.1 

-  6.5  +  jl3 .8 

-109  +  j  0.0 

-1.3  +  J22.1 

-  6.5  -  J13.8 

-10®  +  j  0.0 

-1.3  -  J22.1 

-  6.5  -  j 13 . 8 

-10®  +  j  0.0 

3 .  DISCRETE  DESIGN 

The  approach  takes  into  account  that  the  autopilot  is  implemented 
digitally  in  a  microprocessor.  The  autopilot  update  rate  is  10  ms  and 
there  is  a  6  ms  computational  delay  between  the  autopilot  commands  and 
the  actuator  commands .  The  computational  delay  and  update  rate  are 
shown  in  Figure  11.  In  the  nonlinear  6-DOF  missile  simulation  the 
control  commands  from  the  autopilot  are  delayed  by  6  ms  to  simulate  a 
digital  implementation.  The  continuous  linear  system  is  discretized 
using  the  method  in  Reference  13.  Following  is  a  description  of  the 
method  used  to  discretize  the  continuous  system.  The  control,  u(k) ,  is 
specified  as 


u(k)-  -  Kx(k) 


(117) 


where  u(k)  is  delayed  by  the  computational  delay  r.  In  Figure  11,  the 
transition  from  k  to  k+r,  total  time  equal  to  r,  is  given  by 
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x(k+r)  -  $(r )x(k)  +  r(r)u(k-l) 


(118) 


The  transition  from  k+T  to  k+1,  total  time  equal  to  T-r,  is  given  by 


x(k+l)  -  *(T-r)x(k+r)  +  T(T-r)u(k)  (119) 


Substituting  (118)  for  x(k+r)  into  (119)  leads  to 

x(k+l)  -  $(T-r)*(r)x(k)  +  <&(T-r)r(r)u(k-l)  +  T(T-r)u(k)  (120) 


which  may  be  simplified  to 

x(k+l)  -  $(T)x(k)  +  $(T-r)r(r)u(k*l)  +  T(T-r)u(k)  (121) 

To  remove  u(k-l)  from  the  right  hand  side  of  Equation  (121),  the  state 
vector  xi  is  augmented  with  u(k-l),  this  leads  to  the  discrete  state 
vector 


x 


d 


xx(k) 

u(k-l) 


(122) 


The  discrete  system  in  Equation  (121)  can  now  be  written  as: 


x(k+l) 

$(T)  ®(T-r)T(r) 

(k) 

r(T-r) 

- 

+ 

u(k) 

0  0 

u(k-l) 

I 

u(kT) 


(123) 
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Xr(k+1) 


b - x - *1  I* — x — ►  |  | 

kT  (k  +  1)T  (k  +  2)T 


TIME 

T  =  0.010  Second 

T  =  0.006  Second 
T  -  X  =  0.004  Second 

Figure  1 1 .  Computational  Delay 
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Using  this  approach  the  linear  model  is  augmented  with  the  control 
commands  for  u(k-l).  The  previous  control  commands  are  used  since  these 
are  known,  and  there  are  no  direct  measurements  available  for  the 
controls.  This  is  possible  since  the  control  actuators  response  time  is 
faster  than  the  autopilot  regulator.  This  adds  three  states,  so  the 
discrete  model  now  consists  of  10  states.  During  the  discrete  LQR 
design  no  weighting  is  placed  on  the  augmented  states  for  the  controls, 
the  state  weighting  matrix  Q  is  augmented  with  zeros.  At  this  point, 
the  new  state  vector  is  defined  as 


x-  aflpqrapS  5 
D  L  I  *1  P  q 


(124) 


4.  GAIN  SCHEDULING 

a.  Implementation  Considerations 

An  important  consideration  of  the  design  is  to  minimize  the  amount 
of  gain  scheduling  and  have  large  gain  margins  in  the  design.  Since  the 
model  is  linearized  about  specific  flight  conditions,  gains  must  be 
scheduled  to  account  for  the  other  than  nominal  flight  conditions .  Gain 
scheduling  can  become  an  enormous  burden  on  the  digital  implementation 
of  the  autopilot  so  it  is  desirable  to  limit  the  amount  of  gains  that 
must  be  scheduled. 

A  brute  force  gain  scheduling  technique  is  too  complex  to  implement 

and  includes  too  many  rapidly  changing  variables.  For  the  problem  of  a 

constant  mass  model,  and  a  discretized  model,  the  gains  may  need  to  be 

scheduled  against  the  following  variables:  a,  jj,  p,  q,  r ,  6  ,  S  ,  8  , 

p  r 

mach,  and  altitude.  This  indirectly  includes  dynamic  pressure  since 
dynamic  pressure  is  a  function  of  altitude  and  mach.  There  are  10 
variables  and  3  controls  which  means  that  a  total  of  30  gains  will  have 
to  be  scheduled.  If  mass  were  not  assumed  to  be  constant,  it  too  would 
become  part  of  the  overall  scheduling  problem. 

Normally,  autopilot  gains  are  implemented  by  scheduling  them  as  a 
function  of  only  a  few  variables.  Typically,  these  are  combinations  of 
a,  0,  dynamic  pressure,  and  mach.  While  other  combinations  may  be  used, 
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the  scheduling  is  restricted  still  to  a  small  number  of  variables.  Some 
general  rules  of  thumb  can  be  applied  to  the  gain  scheduling  problem. 

First,  it  is  desirable  to  schedule  gains  against  slower  changing 
parameters  such  as  mach  and  dynamic  pressure.  This  helps  to  maintain 
stability  during  the  presence  of  rapidly  changing  guidance  loop  dynamics 
(primarily  near  intercept) . 

Second,  since  the  implementation  is  discrete,  new  updates  for  gains 
will  be  introduced  as  small  steps  into  the  system  resulting  in  some 
possible  undesirable  transient  responses.  This  means  it  can  be 
necessary  to  slow  the  rate  at  which  the  gains  are  calculated. 

Third,  it  is  possible  to  scale  some  gains  against  a  slowly  changing 
parameter  such  as  dynamic  pressure.  Other  relationships  between  the 
gains  and  variables  if  identified  can  be  exploited.  Gains  that  only 
undergo  small  variations  can  be  fixed  as  constant.  Gains  that  have 
little  contribution  can  possibly  be  ignored. 

Finally,  large  gains  margins  are  desirable  for  the  designs  at  the 
nominal  conditions  so  the  system  is  likely  to  be  stable  when  it  is  not 
at  the  nominal  operating  points.  This  helps  in  terms  of  reducing  the 
amount  of  gains  that  must  be  scheduled.  The  goal  set  for  the  nominal 
design  is  to  have  gain  margins  of  at  least  6dB  and  phase  margins  greater 
than  30  degrees. 


b .  Gain  Schedule . 

The  gain  scheduling  approach  used  here  involves  the  random 
generation  of  500  linear  models  over  a  number  of  flight  conditions, 
using  the  nominal  flight  condition  as  the  mean  of  a  multivariate  uniform 
distribution.  The  mean  components  and  distribution  endpoints  are  shown 
in  Table  3.  Figure  12  summarizes  the  gain  scheduling  approach.  For 
each  randomly  selected  flight  condition  a  discrete  LQR  design  is 
performed  based  on  the  resulting  linear  model.  The  discrete  LQR  designs 
are  then  evaluated  to  assess  the  stability  of  the  resulting  closed  loop 
controller  for  each  design.  Redesign  is  carried  out  if  required.  For 
this  particular  application,  all  designs  were  found  to  be  stable  and  met 
the  desired  gains  and  phase  margins. 
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.  Gain  Scheduling  Approach 


TABLE  3.  FLIGHT  CONDITIONS  DISTRIBUTIONS  USED 
IN  THE  INITIAL  REGULATOR  DESIGN 


VARIABLE 

MEAN 

+ 

a  (degrees) 

10 

15 

$  (degrees) 

0 

8 

p  (degrees/second) 

0 

800 

q  (degrees/second) 

0 

400 

r  (degrees/second) 

0 

400 

altitude  (feet) 

40,000 

0 

mach 

2.5 

0 

5.  MACH  2.5,  AND  40,000  FOOT  DESIGN 

As  shown  in  Table  3,  the  nominal  flight  conditions  are  at  mach  2.5, 
40,000  feet  altitude,  and  angle-of-attack  of  10  degrees.  A  discrete  LQR 
design  for  the  system  linearized  about  the  nominal  flight  condition  was 
performed  using  the  previously  defined  approach  to  specify  the  weighting 
matrices  Q  and  R. 

Figures  13,  14,  and  15  are  discrete  bode  plots  for  the  nominal 
design.  The  frequencies  used  to  generate  the  plots  varied  from  0.1 
radians/second  to  1000  radians/second.  The  plots  are  for  the  augmented 
discrete  open  loop  system  with  the  output  specified  as 


u  -  K  xd  (125) 

Where 

T 

u  -  |?p  8^  (126) 

Figure  13  has  for  its  input  and  output,  5  and  8  .  Figure  14  has  for  its 

—  p  p  ’ 

input  and  output,  8  and  8  .  Figure  15  has  for  its  input  and  output, 
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OS 


Figure  1 4.  Discrete  Bode 


System  1  (Rudder  Loop) 


6  and  S  .  Based  on  these  plots  it  can  be  seen  that  the  phase  and  gain 
margins  for  the  nominal  design  have  been  met. 

The  next  step  was  to  assess  the  designs  over  the  entire  range  of 
conditions  shown  in  Table  3.  Figures  16  through  18  are  plots  for  the 
gains  based  on  the  conditions  in  Table  3 .  Each  plot  was  obtained  by 
plotting  the  gains  against  the  individual  flight  conditions  in  an 
attempt  to  observe  any  relationship  between  gains  and  the  flight 
conditions.  After  a  few  tries  it  readily  became  apparent  that  most 
gains  are  primarily  dependent  on  a  or  p .  A  few  gains  showed  some 
dependence  on  r  but  these  gains  were  small.  In  some  cases  there  is  a 
strong  dependence  of  the  gains  on  both  a  and  p.  Following  is  a  brief 
summary  of  Figures  16  -  18. 

In  Figure  16  (B,  C,  E,  H,  and  J),  Figure  17  (A,  D,  I,  and  J),  and 
Figure  18  (B,  C,  E,  H,  and  J)  a  strong  relationship  between  a  and  the 
gains  appears  to  exist.  Figures  16  G  and  18  G,  show  some  correlation  to 
a.  Very  little  correlation  could  be  discerned  between  any  of  the 
variables  and  Figures  16  (G)  and  18  (G) . 

Figure  16  (F) ,  Figure  17  (F) ,  and  Figure  18  (F)  all  show  a  strong 
relationship  to  the  roll  rate  (p) .  In  Figure  16  (A,  D,and  I),  Figure  17 
(B  and  E) ,  and  Figure  18  (A  and  D)  although  plotted  against  one 
variable,  it  appeared  during  plotting  the  gains,  that  a  relationship 
exists  between  these  gains  and  a  and  p.  In  the  plots  where  a 
relationship  appears  to  exist  between  the  gain  and  the  yaw  rate  (r) , 
Figure  17  (C,  G,  H,  and  J),  and  Figure  18  (D  and  I),  it  also  appears 
that  these  gains  when  multiplied  by  their  corresponding  variables  will 
be  negligible. 

It  is  interesting  to  note  that  the  gains  that  vary  as  a  function  of 
p,  or  a  combination  of  a  and  p,  tend  to  vary  linearly  with  p.  This  is 
not  readily  obvious  from  some  of  the  figures  presented;  however,  while 
plotting  these  gains  against  different  variables  the  relationship  does 
occur.  This  is  a  useful  result  since  p  is  a  rapidly  changing  variable. 
Now  it  appears  that  scheduling  against  p  may  not  pose  problems.  The 
only  real  exception  to  this  is  Figure  17  (F);  however,  in  this  figure 
the  gain  appears  very  close  to  having  a  linear  relationship  with  roll 
rate.  To  isolate  and  identify  these  gains,  200  designs  were  performed 
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using  the  conditions  in  Table  4.  Results  of  these  designs,  although  not 
shown,  verify  the  assumptions  made  about  Figures  16  through  18. 


TABLE  4.  FLIGHT  CONDITIONS  DISTRIBUTIONS  USED 
IN  THE  FINAL  REGULATOR  DESIGN 


VARIABLE  MEAN  ± 

a  (degrees) 

10 

15 

0  (degrees) 

0 

0 

p  (degrees/second) 

0 

800 

q  (degt  os/second) 

0 

0 

r  (degrees/second) 

0 

0 

altitude  (feet) 

40,000 

0 

mach 

2.5 

0 

From  the  designs  based  on  Table  4,  it  is  possible  to  assume  that 
for  a  given  altitude  and  mach  the  gains  need  only  be  scheduled  against  a 
and  p.  In  most  cases  the  gains  need  only  be  scheduled  against  a.  This 
considerably  simplifies  the  approach  to  gain  scheduling.  Now  simple 
table  look-ups  can  be  used  for  gain  scheduling  where  the  independent 
variables  are  a  and  p. 

Another  result  obtained  from  Figures  16  through  18  is  the 
interdependence  of  the  roll  and  yaw  dynamics,  and  independence  of  the 
pitch  dynamics  from  the  yaw  and  roll  dynamics.  This  is  an  interesting 
result  since  autopilot  studies  usually  begin  with  simplifying  the 
nonlinear  model  and  then  linearizing  the  simplified  model.  For  the 
EMRAAT  airframe,  approaches  have  either  separated  the  model  into  pitch 
and  yaw/roll  dynamics  or  separated  the  model  into  roll  and  pitch/yaw 
dynamics.  From  the  figures  it  appears  that  a  reasonable  approach  is  to 
separate  the  model  into  pitch  dynamics  and  yaw/roll  dynamics.  The  model 
for  the  yaw/roll  dynamics  should  possibly  include  a  and  ajas  states. 

The  approach  for  this  problem  does  not  separate  the  dynamics  of  the 
model.  Separating  the  dynamics  of  the  model  is  a  potential  topic  for 
further  work  since  it  leads  to  a  simpler  implementation. 
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Figure  16.  Correlation  Plots  for  Aileron  Gains 


Figure  16.  Correlation  Plots  for  Aileron  Gains  (Continued) 


Figure  1 6.  Correlation  Plots  for  Aileron  Gains  (Concluded) 


Figure  1 7.  Correlation  Plots  for  Elevator  Gains 
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Figure  17.  Correlation  Plots  for  Elevator  Gains  (Continued) 


0.2 
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Figure  17.  Correlation  Plots  for  Elevator  Gains  (Concluded) 


Figure  18.  Correlation  Plots  for  Rudder  Gains  (Continued) 


To  assess  the  importance  of  the  gains  on  a  and  a^for  the  aileron 
and  rudder  commands,  both  a  continuous  and  discrete  design  for  the 
nominal  flight  condition  were  performed  with  the  roll-rate  changed  to 
600  degrees/sec.  First,  the  closed  loop  eigenvalues  for  both  the 
continuous  and  discrete  systems  were  computed.  Next,  the  aileron  and 
elevator  gains  on  a  and  a  were  set  to  zero.  With  the  new  set  of  gains, 
the  closed  loop  eigenvalues  were  computed.  Table  5  contains  these 
results . 

As  can  be  seen  in  Table  5,  the  aileron  and  rudder  gains  on  a  and 
a^appear  to  have  a  significant  impact  on  the  continuous  closed  loop 
eigenvalues.  The  discrete  closed  loop  eigenvalues  do  show  that  one  of 
the  eigenvalues  has  moved  closer  to  the  unit  circle  for  the  case  when 
the  two  gains  were  forced  to  zero.  In  the  gain  schedules  all  of  the 
gains  are  retained  and  none  are  forced  to  zero . 

To  implement  the  gain  schedules,  the  conditions  in  Table  4  were 
varied  using  an  increment  of  1  degree  for  alpha  and  an  increment  of  400 
degrees/sec  for  p.  This  leads  to  155  point  designs.  For  simplicity 
these  designs  are  scheduled  using  simple  2 -dimensional  linear  table 
look-ups . 
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TABLE  5.  COMPARISON  OF  EIGENVALUES 


6. 


MACH  3.0,  AND  65,000  FOOT  DESIGN 

The  autopilot  was  redesigned  using  the  previously  described 
approach,  but  now  for  mach  3.0,  at  65,000  feet.  For  the  nominal  design, 
the  linear  system  in  Figure  9  is  used.  This  leads  to  a  new  state 
weighting  matrix  Q.  A  new  control  weighting  matrix  was  also  specified 
as 


6.2 

2  0 

1.8 

0 

2.5 

0 

CO 

H 

8  0 

8.2 

(127) 


Table  6  contains  the  discrete  open  and  closed  loop  eigenvalues  for  the 
system  of  Equations  (93)  and  (94).  The  same  procedure  for  gain 
scheduling  was  performed  as  was  done  for  the  mach  2.5,  40,000  foot  case. 
When  plotting  gains  little  or  no  change  was  observed  in  the  way  the 
gains  were  correlated  to  a  and  p.  There  is  however,  a  little  more 
dispersion  of  the  points  in  some  of  the  plots.  This  dispersion  can  be 
attributed  to  the  higher  altitude,  which  leads  to  a  lower  dynamic 
pressure.  The  lower  dynamic  pressure  reduces  the  effectiveness  of  the 
controls  and  reduces  the  magnitude  of  the  aerodynamic  derivatives.  Once 
this  happens,  coriolis  and  gyroscopic  coupling  effects  become  more 
pronounced. 
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TABLE  6.  DISCRETE  EIGENVALUES  FOR  SYSTEM  NUMBER  2 


DISCRETE 

SYSTEM  OPEN  LOOP  EIGENVALUES 

0.00 

+ 

jO.OO 

0.00 

+ 

jO.OO 

0.00 

+ 

jO.OO 

0.99 

+ 

jO.  13 

0.99 

- 

jO .  13 

1.00 

+ 

jO .  06 

1.00 

- 

jO .  06 

1.00 

+ 

jO.OO 

1.00 

+ 

jO.OO 

1.00 

+ 

10. OC 

DISCRETE 

SYSTEM  CLOSED  LOOP  EIGENVALUES 

0.00 

+ 

jO.OO 

0.00 

+ 

jO.OO 

0.00 

+ 

jO.OO 

0.81 

+ 

jO.OO 

0.95 

+ 

jO.OO 

0.95 

+ 

jO.08 

0.95 

- 

jO.08 

0.98 

+ 

jO.OO 

0.98 

+ 

J0.06 

0.98 

• 

J0.06 
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SECTION  VII 

TESTING  AND  EVALUATION 


1.  AUTOPILOT  USED  FOR  COMPARISON 

For  purposes  of  comparison,  the  autopilot  in  Reference  7  is  used. 
The  chief  reason  for  this  is  both  the  autopilot  in  Reference  7  and  the 
autopilot  in  this  study  use  the  same  nonlinear  simulation  and  the  same 
Monte-Carlo  simulation.  The  autopilot  in  Reference  7  will  be  referred 
to  as  the  baseline  autopilot. 


2 .  ACCELERATION  STEPS 

Step  commands  for  specified  accelerations  in  the  inertial  reference 
frame  were  performed  as  shown  in  Table  7  for  the  mach  2.5,  40,000  foot 
design.  This  approach  is  a  good  method  to  assess  autopilot  performance 
since  an  inertial  guidance  step  demands  more  from  the  autopilot  than  a 
step  on  one  of  the  states  in  x.  It  should  be  noted  that  more  stressing 
maneuvers  should  be  used  for  a  complete  analysis  of  performance,  but  for 
this  study,  only  the  four  steps  shown  in  Table  7  are  used.  In  all  four 
cases  the  missile  initially  assumes  straight  and  level  flight  with  an 
angle-of-attack  that  is  sufficient  to  maintain  trim  (keep  altitude 
constant).  The  steps  are  initiated  at  1  second  into  flight. 

In  step  number  one,  the  vehicle  is  commanded  to  accelerate  in  only 
the  inertial  y-plane.  Step  1  is  the  least  complex  maneuver  of  the  four 
steps,  however,  the  missile  must  bank  about  its  velocity  vector.  To  do 
this,  the  desired  a  must  be  achieved,  and  ft  must  be  kept  to  a  minimum. 

In  step  2,  the  missile  is  commanded  to  perform  a  10  g  dive.  This  forces 
the  missile  to  bank  180  degrees  about  its  velocity  vector.  Step  3 
combines  steps  1  and  2.  In  step  4,  the  missile  is  commanded  to 
accelerate  in  the  inertial  y  and  z  planes.  This  maneuver  causes  the 
missile  to  climb.  Note  that  the  step  response  will  be  less  than  the 
step  command,  this  is  a  result  of  the  BTT  logic.  Also,  it  should  be 
noted  that  the  acceleration  step  responses  will  lose  magnitude  in  steady 
state,  this  is  due  to  the  fact  that  the  airframe  cannot  accelerate 
indefinitely . 
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TABLE  7.  INERTIAL  GUIDANCE  STEPS 
(MACH  2.5,  40,000  FT  ALTITUDE) 


STEP  # 

STEPS  (g's) 

a  a 

y  z 

J  c  c 

1 

25 

0 

2 

0 

25 

3 

25 

25 

4 

25 

-25 

Figure  19  shows  the  response  to  step  1.  The  solid  line  in  Figure 

19  is  the  achieved  acceleration  for  a  and  the  dashed  line  is  the 

yi 

achieved  acceleration  for  a  .  From  Figure  19 ,  the  time  response  for 

*i 

the  step  1  command  to  67  percent  of  the  steady  state  value  is 

approximately  0.5  seconds.  This  is  considerably  slower  than  the 

baseline  autopilot  which  has  a  response  of  0.3  seconds.  The  transient 

response  for  a  is  less  than  that  of  the  baseline  autopilot.  In  Figure 
*i 

19  there  appears  to  be  some  instability  in  the  a  response  at  1.1 

yi 

seconds.  This  is  not  present  in  the  response  for  the  baseline 

autopilot.  The  a  response  shows  a  bias  slightly  less  than  1.  This 
*i 

bias  is  necessary  to  offset  gravity. 

Figure  20  shows  the  step  2  response.  The  autopilot  achieves  67 
percent  of  steady  state  at  0.7  seconds.  The  baseline  autopilot  achieved 
steady  state  at  0.6  seconds.  The  a  transient  response  in  the  baseline 

yi 

autopilot  is  much  more  severe,  with  a  peak-to-peak  magnitude  of  8  g's. 

Figure  21  shows  the  step  3  response.  The  autopilot  achieves  67 

percent  of  steady  state  in  0.6  seconds  for  both  a  and  a  The 

yi  *i 

baseline  acceleration  responses  were  0.4  seconds  and  0.5  seconds 
respectively.  The  baseline  response  has  more  transients  at  1.1  seconds 
and  after  1.4  seconds. 
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Figure  22  shows  the  step  4  response.  Both  inertial  accelerations 
achieve  67  percent  of  steady  state  in  0.3  seconds  and  the  baseline 
autopilot  response  takes  0.2  seconds.  Figures  23,  24,  25,  26,  and  27 
show  a,  fi,  p,  q,  and  r  responses  for  step  4.  The  response  for  a  is  0.2 
seconds  less  than  the  baseline  a  response.  There  are  some  other  minor 
differences  in  the  other  responses  shown  and  those  of  the  baseline 
autopilot;  however,  it  appears  that  the  a  response  difference  is  the 
most  critical  since  it  is  significantly  slower  than  the  baseline. 

In  addition  to  the  steps  shown  in  Table  7,  100  g  steps  were  also 
used  to  assess  acceleration  capability  of  the  missile.  Based  on  these 
results,  the  maximum  acceleration  of  the  airframe  is  slightly  less  than 
30  g's. 


3.  MISSILE -TARGET  INTERCEPTS 

Three  engagement  conditions  were  used  to  test  the  autopilot 
designs.  These  engagement  conditions  were  used  for  both  the 
deterministic  runs  and  the  Monte  Carlo  runs.  Table  8  shows  the  missile 
engagement  conditions  used,  and  Table  9  shows  the  target  engagement 
conditions  used. 


TABLE  8.  MISSILE  ENGAGEMENT  CONDITIONS 


MISSILE 

ENGAGEMENT 

ALTITUDE 

MACH 

1 

40,000  ft 

2.5 

2 

40,000  ft 

2.5 

3 

65,000  ft 

3.0 
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TABLE  9.  TARGET  ENGAGEMENT  CONDITIONS 


TARGET 

ENGAGEMENT 

ALTITUDE 

MACH 

RANGE 

HEADING 
(degrees  off 
head-on) 

1 

40,000  ft 

2.0 

15,000  ft 

15 

2 

40,000  ft 

1.1 

15,000  ft 

180 

3 

65,000  ft 

2.3 

15,000  ft 

5 

4.  MONTE  CARLO  RUNS 

Miss  distance  sensitivity  of  the  autopilot  to  modeling  errors  and 
measurement  noise  was  investigated  to  determine  the  effects  these  have 
on  the  miss  distances.  For  the  purpose  of  comparison,  the  Monte  Carlo 
runs  emulated  the  runs  accomplished  for  the  autopilot  design  of 
Reference  7.  The  number  of  runs  per  Monte  Carlo  simulation  was  100. 
Numerous  error  sources  are  included  in  the  Monte  Carlo  simulation. 
Included  in  these  were  sensor  noise,  sensor  scale  factor  and  bias 
errors,  variations  on  the  deviation  from  the  initial  conditions,  and 
variations  on  wind  angle  and  velocity.  Also  included  were  variations  on 
the  aerodynamic  parameters  to  simulate  modeling  errors. 

In  the  Monte  Carlo  simulation  the  noise  sequences  used  for  the 
measurements  are  characterized  by  a  normal  distribution.  At  the 
beginning  of  each  run  a  random  draw  was  performed  to  establish  the  mean 
and  standard  deviation  for  each  of  the  noise  sequences  used.  The  random 
draws  themselves  are  each  characterized  by  a  normal  distribution,  and 
are  established  as  an  input  to  the  simulation. 

Tables  10  and  11  contain  a  brief  summary  of  the  Monte  Carlo  results 
for  engagements  1  and  2.  The  assumption  used  in  constructing  Tables  1 
and  2  is  that  100  runs  per  Monte  Carlo  simulation  were  enough  to  treat 
the  miss  distances  as  if  they  were  normal  distributions.  Only  miss 
distances  of  50  feet  or  less  were  considered.  Performance  in  engagement 
1  is  only  slightly  better  than  the  performance  of  the  baseline 
autopilot.  Performance  in  engagement  2  is  better  than  the  baseline 
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autopilot  performance.  Performance  of  the  65,000  feet,  mach  3.0  design 
in  engagement  3  was  not  compared  to  the  baseline  since  the  Monte  Carlo 
results  necessary  for  a  direct  comparison  were  unavailable.  In 
engagement  3,  the  mean  miss  distance  was  22.3,  with  a  standard  deviation 
of  10.3. 


TABLE  10.  MONTE  CARLO  RESULTS  FOR  ENGAGEMENT  1 


MISS  DISTANCE 

SUMMARY 

Mean 

Standard 

Deviation 

40,000,  mach  2.5  design 

9.74 

8.24 

Baseline  autopilot 

10.40 

8.22 

TABLE  11.  MONTE  CARLO  RESULTS  FOR  ENGAGEMENT  2 


MISS  DISTANCE 

SUMMARY 

Mean 

Standard 

Deviation 

40,000,  mach  2.5  design 

12.10 

11.51 

Baseline  autopilot 

16.50 

11.63 
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Figure  19.  Response  for  Step  1 
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Figure  20.  Response  for  Step  2 
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Figure  21 .  Response  for  Step  3 
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Figure  22.  Response  for  Step  4 
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igure  23.  Angle  of  Attack  Response  for  Step  4 
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Figure  24.  Sideslip  Response  for  Step  4 
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Figure  25.  Roll  Rate  Response  for  Step  4 
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Figure  26.  Pitch  Rate  Response  for  Step  4 
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Figure  27.  Yaw  Rate  Response  for  Step  4 


SECTION  VIII 

CONCLUSIONS  AND  RECOMMENDATIONS 


1.  CONCLUSIONS 

An  autopilot  was  designed  based  on  the  full  nonlinear  missile 
dynamics.  Performance  of  the  autopilot  was  evaluated  with  Monte  Carlo 
simulations.  Although  deterministic  performance  of  the  autopilot  is 
slightly  worse  than  performance  of  the  baseline  autopilot,  the 
Monte  Carlo  simulation  of  the  autopilot  demonstrated  improved 
performance  over  the  baseline  autopilot.  The  autopilot  response  is  more 
damped  than  the  baseline  autopilot  response.  The  baseline  autopilot 
leads  to  more  coupling  in  the  acceleration  step  responses. 

The  improved  Monte  Carlo  results,  and  the  slower  step  responses  of 
the  autopilot  as  compared  to  the  baseline  autopilot,  are  due  to  the 
tradeoff  done  for  selection  of  the  control  weighting  matrix  (R)  used  in 
the  regulator  design.  An  improved  step  response  leads  to  improved 
deterministic  performance,  but  the  autopilot  becomes  less  robust  in  the 
presence  of  sensor  noise  and  more  sensitive  to  modeling  errors.  For 
these  reasons  R  was  selected  to  improve  the  Monte  Carlo  results. 

The  approach  to  the  design  also  demonstrated  a  way  to  assess  model 
cross  coupling.  For  the  EMRAAT  airframe,  this  approach  demonstrated 
that  for  decoupled  approaches,  decoupling  the  pitch  dynamics  from  the 
roll  and  yaw  dynamics  appears  to  be  good  approach.  The  approach 
demonstrated  the  coupling  between  the  roll  and  yaw  dynamics.  The 
approach  also  demonstrated  that  some  pitch  coupling  is  present  in  the 
roll  and  yaw  dynamics.  The  65,000  foot,  mach  3.0  design,  showed  only 
slightly  more  coriolis  and  gyroscopic  coupling  than  the  designs  at 
40,000  foot,  mach  2.5.  In  either  case,  the  coriolis  and  gyroscopic 
coupling  did  not  pose  design  or  implementation  problems. 

The  stochastic  approach  used  to  assess  gain  scheduling  requirements 
demonstrated  that  the  gain  scheduling  requirements  for  the  autopilot  are 
minimal  and  may  be  scheduled  as  a  function  of  angle -of -attack  and 
roll  rate.  Leaving  roll  rate  out  of  gain  scheduling  degrades 
performance,  but  this  is  a  possible  tradeoff  if  simpler  gain  schedules 
are  required.  The  exceptions  to  this  are  the  gains  that  vary  as  a 


87 


function  of  roll  rate  and  couple  pitch  dynamics  into  the  yaw  and  roll 
controls,  these  gains  should  not  be  neglected. 


2 .  RECOMMENDATIONS 

a.  Model  Order  Reduction 

Based  on  the  results  of  the  approach  used  in  this  study,  a  next 
generation  autopilot  design  might  center  around  a  pitch  model  and  a 
roll/yaw  model.  The  roll/yaw  model  would  include  the  coupling  from  the 
pitch  channel.  The  advantage  to  this  approach  is  it  results  in  a 
reduced  set  of  gains  to  schedule.  This  will  also  allow  comparison  of 
this  model  to  the  full  order  model  to  identify  the  real  benefits  of 
either  approach. 

b.  More  Dynamic  Flight  Envelope 

The  approach  should  be  extended  to  a  varying  flight  envelope  with 
larger  altitude  and  mach  variations,  which  leads  to  larger  dynamic 
pressure  variations.  A  thrusting  motor  should  be  included  to  add  the 
effects  of  changing  missile  mass. 


c.  Reducing  Sensitivity  to  Noise  and  Missile  Flexible  Modes 
A  major  drawback  of  the  approach  used  in  this  study,  was  the 
inability  to  improve  both  deterministic  performance  and  performance  in 
the  Monte-Carlo  simulations.  Since  the  autopilot  is  working  in  an 
environment  of  noisy  measurements  and  modeling  errors,  techniques  such 
as  LQG/LTR  should  be  employed  to  reduce  sensitivity  to  noise.  Analysis 
should  be  extended  to  include  body  flexible  modes.  Notch  filters,  for 
example ,  may  be  employed  to  attenuate  bending  modes . 
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THE 

APPENDIX 

EMRAAT  AERODYNAMIC  DATA 


89/90  (Blank) 


This  appendix  lists  the  data  used  for  the  extended  medium  range 
air-to-air  technology  (EMRAAT)  airframe.  The  data  listed  is  for  the 
airframe  with  no  fuel  (empty  weight).  The  first  variables  listed  are 
those  that  do  not  change. 


d  Missile  Reference  Diameter  (0.625  ft) 
S  Missile  reference  area  (0.3067  ft2) 

W  Missile  Weight  (empty)  (227  lb) 


Moments  of  Inertia  (slug  -  ft  )  : 


I  -  1.08 
xx 

I  -  70.13 

YY 


zz  -  70.66 


Products  of  Inertia  (slug  -  ft  ) 


xy  -  0.274 

Xxz  -  0 . 704 

^z  -  -0.017 

The  following  data  is  a  function 

of  the  flight 

data  is  listed  for 

two  conditions: 

Condition  #1 : 

Angle-of-attack 

10  degrees 

Sideslip  angle 

0  degrees 

Roll  rate 

0  deg/sec 

Pitch  rate 

0  deg/sec 

Yaw  rate 

0  deg/sec 

Altitude 

40,000  ft 

Mach 

2.5 

Condition  #2: 

Angle-of-attack 

10  degrees 

Sideslip  angle 

0  degrees 

Roll  rate 

0  deg/sec 

Pitch  rate 

0  deg/sec 

Yaw  rate 

0  deg/sec 

Altitude 

65,000  ft 

Mach 

3.0 

The 
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Q  - 


2 

Dynamic  Pressure  (lb/ft  )  : 

Condition 

Condition 


1: 

2: 


1719.93 

749.29 


Following  are  the  aerodynamic  derivatives  (units  are  per  rad) . 

Note  that  C  was  obtained  by  taking  the  slope  of  the  aerodynamic 
R 

a 

derivative  C^  (not  listed) . 


Condition  1 

Condition  2 

c 

N 

a 

-50.7067 

-47.3550 

C 

a 

0.0221 

0.0182 

C 

N 

q 

-0.0103 

-0.0082 

C 

"«q 

-3.1182 

-2.4653 
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Condition  1 


Condition  2 


-16.1902 

-16.2871 

-0.0004 

-0.0003 

0.0082 

0.0051 

0.2229 

0.5171 

3.5252 

3.2733 

Condition  1 


Condition  2 


-9.3923 

-7.5674 

-0.0069 

-0.0027 

0.0011 

0.0051 

-4.0739 

-3.2080 

-5.5816 

-4.7366 
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Condition  1 

Condition  2 

c 

m 

a 

-17.6108 

-16.8330 

n 

u 

m* 

a 

0.0124 

0.0182 

c 

m 

-0.1329 

-0.1043 

C 

-25.4394 

-20.4340 
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